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Abstract 

In order to have massive neutrinos, the right-handed neutrino/sneutrino superfield (N) need to 
be introduced in supersymmetry. In the framework of NMSSM (the MSSM with a singlet S) such 
an extension will dynamically lead to a TeV-scale Majorana mass for the right-handed neutrino 
through the SNN coupling when S develops a vev (the free Majorana mass term is forbidden by the 
assumed Z3 symmetry). Also, through the couplings SNN and SH,H,, the SM-like Higgs boson 
(a mixture of H,, Hg and S) can naturally couple with the right-handed neutrino/sneutrino. As 
a result, the TeV-scale right-handed neutrino/sneutrino may significantly contribute to the Higgs 
boson mass. Through an explicit calculation, we find that the Higgs boson mass can indeed be 
sizably altered by the right-handed neutrino/sneutrino. Such new contribution can help to push 


up the SM-like Higgs boson mass and thus make the NMSSM more natural. 


PACS numbers: 12.60.Jv,11.30.Qc,12.60.Fr,14.80.Cp 


I. INTRODUCTION 


Supersymmetry (SUSY) |1, 2] gives a natural solution to the hierarchy problem suffered 
by the Standard Model (SM). Also, it provides a good dark matter candidate and realizes 
the gauge coupling unification. Among the SUSY models, the Minimal Supersymmetric 
Standard Model (MSSM) [3] has been intensively studied. However, the recently discovered 
Higgs-like boson around 125 GeV caused a problem for this model, i.e., a 125 GeV Higgs 
boson requires a heavy stop or a large tri-linear coupling A; and thus incurs the little 
hierarchy problem. Besides, the MSSM suffers from the -problem [4]. 

It is remarkable that both the little hierarchy problem and the p-problem can be solved in 
the Next-to- Minimal Supersymmetric Standard Model (NMSSM) [5] in which an additional 
gauge singlet S is introduced (in fact, the NMSSM was proposed even earlier than the 
MSSM [6]). In this model the j-problem is solved by the dynamical generation of the u- 
term through the coupling S H,H, when S develops a vev, while the little hierarchy problem 
is solved by the generation of an extra tree-level mass term for the SM-like Higgs boson (thus 
the stop mass or A; is no longer required to be unnaturally large). 

Note that in order to have massive neutrinos, right-handed neutrino/sneutrino super- 
field(s) (N) need to be introduced in SUSY models. For the NMSSM with such right-handed 
neutrino/sneutrino field(s) [7], some intriguing features are present. Due to the assumed Z3 
symmetry, the free Majorana mass term for the right-handed neutrino is forbidden in the 
superpotential. Instead, a TeV-scale Majorana mass for the right-handed neutrino is dy- 
namically generated through the SNN coupling when S develops a vev (note that such 
a TeV-scale majorana mass is too low for the see-saw mechanism and thus the neutrino 
Yukawa couplings H,,LN must be very small). In the same way, a TeV-scale mass for the 
right-handed sneutrino can also be generated, which can serve as a good dark matter can- 
didate [9]. Further, through the couplings SNN and SH,H4, the SM-like Higgs boson (a 
mixture of H,, H4 and S) can naturally couple with the right-handed neutrino/sneutrino. 
As a result, the TeV-scale right-handed neutrino/sneutrino may significantly contribute to 
the Higgs boson mass (in the MSSM with split SUSY, the right-handed neutrino/sneutrino 
can also make sizable contribution to the Higgs mass, as studied in |10]). In this paper we 
will perform an explicit calculation for such contribution. 


This work is organized as follows. In Sec. II we present the spectrum and couplings 


for the Higgs boson and right-handed neutrino/sneutrino. In Sec. III the renormalization 
scheme is described. Numerical results and discussions are given in Sec. IV. Finally, we 


give a summary in Sec.V. 


II. HIGGS AND RIGHT-HANDED NEUTRINO/SNEUTRINO IN NMSSM 
A. Model description 


The NMSSM with a right-handed neutrino superfield N has a superpotential given by 


W = Wynmssm + ANSNN + yn fl, : LN, 
1 
Wyussm = Ya Hu : Qur — YaHa: Qdg — Y.Ha* Leg + ASH,- Hat 359^ (1) 


where the flavor indices are omitted and the dot denotes the SU (2) j, antisymmetric product. 
Since a global Z3 symmetry is imposed, there are no supersymmetric mass terms (like H, Ha, 
NN or SS) in the superpotential. Note that in this model we impose R-parity and thus the 
terms NNN and SSN are forbidden. As a result, the sneutrino-Higgs mixing is avoided 
and also there is no vev for the right-handed sneutrino (we will show how to get the globle 
minimum in the following). Although a bare Majorana mass term N.N is forbidden in the 
superpotential, a TeV-scale Majorana mass can be generated through the coupling SNN 
when S develops a non-zero vev (vs). Such a TeV-scale Majorana mass is too small for the 
conventional see-saw mechanism and thus the Yukawa coupling y. H,L.N should be very 
small (yy <1). Note that here we introduce only one right-handed neutrino superfield to 
illustrate its effects on the Higgs mass. In order to explain the neutrino masses and mixings, 
more right-handed neutrino superfields need to be introduced, each of which will contribute 
to the Higgs mass. In this case, the calculation method is same as in our calculation, but 
the total effects may be more sizable due to more free parameters. 

The soft SUSY breaking terms for Higgs and right-handed sneutrino are given by (here- 


after we use N and N to denote respectively right-handed neutrino and sneutrino) 


-Lot = M}, |Ha]? + MZ |Hal? + M2|SP + (AA, - HaS + SAnS + h.c.) 


-M2|NP + (ANANS NN + h.c.) (2) 


Here we neglected the mixing between left-handed and right-handed sneutrinos because the 


mixing is assumed to be suppressed by yy. In the following we briefly discuss the neutral 


Higgs neutrino sectors. 


B. The neutral Higgs sector 
From Eq. (1) and Eq. (2) we get the Higgs potential 


V = X(|H4*|SP + |HaP|SP + Hu - Hal?) + &^|S?P 
1 
FAR, - HaS*S* + h.c.) +3 g (Hu — |Hal?)? 


1 
-5 gal Hr (Ha)* + HEUS + Mg, HU + Mg, Hal + Ms|S/° 


1 
+AA Ay, - H48 + 3 KA, SÈ + hie.) (3) 


where g? = (g? + g2)/2 with gı and go being the SM gauge coupling constants. Assuming 
H,, H4 and S get vevs such that 
Bs Re(H9)) + ilm(H9) Re( H9) + ilm(H9) Sg 4 i$, 
AL uU /2 9 /2 ? /2 


we can get the mass terms for the Higgs fields, which are presented in [11]. Here we only 


HÌ = v4 + S=ust+ (4) 


show the conventions and give some brief comments: 


1. The mass matrix for the CP-even neutral Higgs is obtained from the real components of 
the Higgs fields. In the basis h*”* = [Re( H9), Re(H9), Sp] and using the minimization 
equations to eliminate the soft masses, one obtains three CP-even states (ordered in 
mass) 


h; — Sihi" (5) 
with an orthogonal rotation 5;;. 


2. The mass matrix for the CP-odd neutral Higgs is obtained form the imaginary com- 
ponents of the Higgs fields [Im(H°), Im(H$), Sz]. Its diagonalization is performed in 
two steps. First, one rotates it into a basis (A, Sr, G) where G = — sin GIm(H°) + 
cos 3Im(H9) is a massless Goldstone mode (tan 3 = v,/va is the ratio of the vevs of 
the two Higgs doublets). Dropping the Goldstone mode, the remaining 2 x 2 mass 


matrix M? in the basis (A, Sr) can be diagonalized by an orthogonal 2 x 2 matrix Pj 


into two physical CP-odd states a; (ordered in mass): 


ay = PA + P51, 


ag = PaA + Po2S7. (6) 


3. The neutralino mass matrix My in the basis V? = (—iM, —ià2, VO, Y9, Ys) can be 
diagonalized by an rotation matrix N;;. Then one obtains five eigenstates (ordered in 


mass) x? = Nijv?. 


C. Right-handed neutrino/sneutrino sector 


Since there is no Dirac mass term here, the mass spectrum of the right-handed neutrino 
sector is very simple. Denoting Ñ = R 4- iM, there are only one CP-even right-handed 
sneutrino (denoted as R) and one CP-odd right-handed sneutrino (denoted as M). The 
right-handed neutrino is denoted as N. From Eq. (1) and Eq. (2) we can get the spectrum 


as 


m? —O4ASUL + Mz + 2r0yusAn + 2An (kv? — Wava) 


me, = AM? + M% — 2\yusAn — 2Àn (KU? — Avava) 


My = 2A NVs. (7) 


With the above spectrum we can get the couplings between the Higgs and the right-handed 
neutrino/sneutrino. Note that in our numerical study we require Mz and Mj, be positive, 
and, as a result, the global minimum of the scalar potential locates at the zero point of 


the right-handed sneutrino field (the right-handed sneutrino has no vev and thus R-parity is 


preserved). In the following we list the couplings which will be used in our later calculations: 


Vig = V2ANA (vuSja + va851) — V2 (2A v, + AM, + AN An) S33, (8) 
Varo = —V2ANA (vuSjs + vaSj) + V2 (2AwRv, — 4ANUs + AN An) S33, (9) 
Vin;RR = Àn 258383 — A(Sj1Si2 + Sa Sj2)] — AM S355, (10) 
Vin;uM = An [259595 — A(Sj1Si2 + 93553)] — 4A 9353, (11) 
Varm = —2An(—Av cos 28P4/V/2 + V2kv, Pi) + V2AN An Pio, (12) 
Vaa;nR = 2Ay(Asin B cos 8P4 Pj + KPi2Pj2) — ANS Pi Pia, (13) 
Va,aj;MM = —2An(Asin 8 cos bPa Pj + &Pi3Pj3) — ANS Pi Pia, (14) 
Van = —V2ANn Sis Vann = V2iAN Pay, (15) 
Vy;RN = —ÀN Lm VMN = E. (16) 


III. RENORMALIZATION SCHEME 


To calculate the neutrino/sneutrino contribution to the Higgs mass, we must calculate 
the one-loop Higgs propagator and choose a renormalization scheme. Here we follow |12] and 
choose the mixed renormalization scheme (other schemes give similar results). We choose 


the following parameter set 


Mz, My, Mus, e, tg, tg, tu, tan, A, Us, K, Ax, (17) 

ORIS Rae DR scheme 
where £j, tg, tg, are the tadpoles of the CP-even Higgs fields. Since we concentrate on 
the right-handed nuetrino/sneutrino contributions, the input parameters from the gauge 
interaction part need not be renormalized. For the parameters which need renormalization, 


we replace them by the renormalized ones plus the corresponding counterterms: 


ty, Oty, -ótg,, tan — tan --ótanf 
tH `> te, de À—A-ÓA 


ty, > tH, + ÔtH,, KoOK+OK 


Us — Us + ÓV,, A, >A, + 6A, . 


In the following we will show how to determine the counter terms in the mixed renormal- 


ization scheme. 
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FIG. 1: Feynman diagrams for the two-point renormalized Higgs functions. 
First, the Higgs doublet and singlet fields are replaced by the renormalized ones: 
Hi >V Zm a= ( + 502.) Hu 
Ha > /Zg, Ha = ( 4 jn) Ha (19) 
So VZs5S = (1+ 502s) s. 


Then the renormalized two-point functions can be obtained from the Feynman diagrams 


shown in Fig. 1 


SH; (k?) = Sik Sj X20) (i,j, k,l = 1, 2,3), (20) 


Maa) eJ 45512) (21) 


where 5;; and Pj; are the matrix elements defined in Eqs.(5) and (6). The renormalization 
condition can be set as 


OY p, gu, (K?) UN 
Ok? k2=(M(?))2 


a 


where M E denotes the corresponding tree-level Higgs mass, and ’div’ shows that we chose 
the DR renormalization scheme which means that in the field renormalization only the 
divergent part A = 2/(4 — D) — yg + In(4v) (yg is the Euler constant) is kept. The field 


renormalization constants ÔZ g}, ÔZ H4, ÔZs are obtained by solving the equations 


óZu,u,; = Sal Za, + |S Zu, + |Si| iZ (¢=1,2,3); (23) 


We use the field renormalization constants to determine the conterterms listed in Eq. (18). 
The detailed calculations are lengthy. In the following we only present the final results and 


give some necessary comments. 


1. Tadpole parameters: 
The tadpole parameters are determined by the condition that they vanish after the 
renormalization. The Feynman diagrams are shown in Fig. 2 and the counter terms 


are determined by 
peu (mde $1529). (24) 


where a" denote the one-loop Higgs tadpoles. 
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FIG. 2: Feynman diagram for the Higgs tadpoles. 
2. The parameter tan f: 
t 
ótanf = | i De cs 22] (25) 
div 
3. The coupling A: 
E È (M2 ) 
i AME, 82, PENTR RIL | ae 


The self-energy 3p; is obtained from the self-energies in the mass eigenstate basis 


DA; A; (i, j = 1,2,3) through 


Spi = Py daa; Pj o. (26) 


4. The singlet Higgs vev v;: 


ÔU; = —Us— ; (27) 


5. The coupling &: 
kK is renormalized through the neutralino renormalization whose diagrams are shown 
in Fig. 3. Note that we have different conventions of vev and thus the formula is a 


little different from Ref. [12]. 


1 OUs 
OK = 2,, (Mw) ss — K (28) 
Uv " 
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FIG. 3: Feynman diagrams for the renormalized two-point neutralino functions. 


6. Tri-linear coupling A,: 


A, is renormalized by the CP-odd Higgs element M, and is given by 


1 
— [Eres (M$3s) - 97] - Axl 


S 


ÔK OU, 
— + 


K Us 


5A, = [- Jos (29) 


where the fuction f can be found in Ref. [12]. 


After the determination of the counterterms, we put these terms into the Higgs mass matrix 
which is shown in the Appendix. Also, by adding the loop contribution to the Higgs mass 


matrix, we can get the mass correction for the Higgs boson. 


IV. NUMERICAL RESULTS 
A. "The right-handed neutrino/sneutrino correction to the Higgs boson mass 


In our numerical calculation we concentrate on the SM-like Higgs boson which is the 


lightest CP-even Higgs boson dominated by the Higgs doublets. From the superpotential in 


FIG. 4: The right-handed neutrino/sneutrino contribution to the SM-like Higgs boson mass versus 
AAN- 


Eq. (1) we can see that the right-handed neutrino/sneutrino interacts with the doublet only 
through the F-term of the singlet Higgs S, and thus the parameter A will play an important 
role in the correction to the Higgs boson mass. Also, from the superpotential we can also see 
that the right-handed neutrino/sneutrino couples to the Higgs sector through the parameter 
Ay. So, as Ay approaching zero, the right-handed neutrino/sneutrino should decouple form 


the NMSSM sector. To check this numerically, we scan the parameter space in the range: 


0< à, k, An «1, 2< tan < 50, (30) 
0< u,Mę <1 TeV, —1 TeV < Ay, Ak, An < 1 TeV, (31) 


Note that in the calculation of the Higgs mass spectrum we chose to use u (= Av,) as an 
input parameter because it is commonly used in the NMSSM phenomenology studies and 
the relevant numerical packages. Also we note that À and Ay may be rather constrained 
(e.g., A at weak scale must be below 0.7) if we require perturbativity of the theory up to 
the grand unification scale |13]. Of course, if we just treat NMSSM as a low energy effective 
theory, such a stringent perturbativity constraint will be much relaxed. 

The correction to the Higgs boson mass versus product of A and Ay is shown in Fig.4. 
From the figure we can see that when the product of A and Ay approaches to zero, the 


correction will approach zero; when AAy is at order 1, the right-handed neutrino/sneutrino 


10 


will alter the mass significantly. Thus, if A and Ay is not small, then the right-handed 
neutrino/sneutrino contribution to the Higgs boson mass must be taken into account. 

Now we check the SUSY limit in the right-handed neutrino/sneutrino sector. From Eq. 
(7) we can see that with Mx and Ay approaching zero, the right-handed neutrino/sneutrino 
sector has a SUSY limit for kv? = Av,vg. In our second scan, we assume the relation 
ku? = Av,vq and let the parameter A, k, tan, Ax, A,, Mg and Ay vary randomly in 
range as in Eqs. (30, 31), only fixing Ay = 0.9. The results are shown in Fig. 5. The results 
show that with , | M3 + Aĉ, approaching zero, the Higgs mass correction approaches zero, 
which confirms the SUSY limit. 


1 E 


0.9 - 4 
E j J 
L X » 4 
d E Aw = 0.9 M E 
E "EX J 
ac 0.7 - ": x mi 
& 0.6 - E 4 
- | V E 
A 0.5 F KM a 
E 04 - a E 
0.3 - 4 
02 - : E 
E go *x 
0.1 - H "m m 
E P Sd 
E x Xx od 
1 


^ 40 10° 
2,A2 
VMZ«AZ (GEV) 
FIG. 5: The right-handed neutrino/sneutrino contribution to the SM-like Higgs boson mass versus 


y Mg + AN. 


It is well known that the Higgs mass can be enhanced by the hierarchy between the 5M 
particles and their SUSY partners. In the right-handed neutrino/sneutrino sector, the mass 
hierarchy between sneutrino and neutrino is controlled by the soft parameters Mx and Ay. 
In order to show the dependence on the mass splitting, we chose a benchmark point: 

Àz0.2, Ay = 0.35, & = 0.4, tang = 10, 
u = 200 GeV, Aj = 300 GeV, A, = —500 GeV, (32) 
and scan the other two parameters in the range of 0 < Mẹ < 1 TeV and —1 TeV < Ay < 


1 TeV. The results are shown in Fig. 6, in which the left panel shows dm; versus Mx and 
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FIG. 6: The right-handed neutrino/sneutrino contribution to the SM-like Higgs boson mass versus 


the sneutrino soft mass Mx, and the ratio m2,/m3, (for mg and my, see Eq. (7)). 


the right panel shows ôm, versus m2,/m3,. From this figure we can see that as Mx increases 
(the mass slitting between sneutrino and neutrino also increases as shown in Eq. (7)), the 
mass correction increases. 

From the above results we see that the right-handed neutrino/sneutrino can either en- 
hance or reduce the Higgs boson mass. Since the parameter space is multi-dimensional (9 
input parameters), we perform an intensive scan to try to figure out what parameter(s) 
determine the sign of the correction. We scan the parameter set (A, Aw, Mg, An) while 
fix other parameters as listed in Eq. (32). The results are shown in Fig. 7. We see that the 
parameter A plays the most important role in this aspect although it cannot solely determine 
the sign. For a large value of A the sign of the correction tends to be negative. Clearly, the 


sign is not sensitive to Ay. We also checked that the sign is not sensitive to other parameters. 


B. Higgs mass with all loop corrections under current experimental constraints 


In the preceding section we only considered the loop corrections from the right-handed 
neutrino/sneutrino. Of course, the loop corrections from other particles (especially the top 
and stop) should also be taken into account. In our following numerical study, we include all 
available loop corrections by using the package NMSSM Tools [11]. Since the right-handed 


neutrino/sneutrino is a gauge singlet, it will not change the Higgs decay or the annihilation 
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0.5 


FIG. 7: The right-handed neutrino/sneutrino contribution to the SM-like Higgs boson mass shown 
in the plane of Ay versus A. Here we scan the parameter set (A, Aw, Mx, An) while fix other 
parameters as listed in Eq. (32). The red ’x’ are for bmp, < —1 GeV, the green ’+’ for -1 GeV < 


dmp < 0, the blue ’o’ for 0 < ômp < 1 GeV, and the magenta ’A’ for dm; > 1 GeV. 


of the dark matter. So, we just add the right-handed neutrino/sneutrino correction to the 
Higgs boson mass in the NMSSM Tools. Then we scan the NMSSM parameter space in the 


range: 


Ü« A, X«L 2< tanß < 50, 
0 < (u, Mı = M2/2 = M3/6, mg, m; = m; = Mz = mg) < 1 TeV, 
—] TeV < (Ax, Ax, At = Ay = A, = Ap) < 1 TeV. (33) 


For the neutrino/sneutrino sector, we set Ay = 0.5 and scan Mx, Aw in the range 
0« M& «1TeV, —1TeV < Ay « 1 TeV. (34) 


In our scan we consider the following experimental constraints [14]: (1) We require the 
lightest neutralino X? to account for the dark matter relic density 0.105 < Qh? < 0.119; 
(2) We require the SUSY contribution to explain the deviation of the muon ay, i.e., až? — 
a»M = (25.5 + 8.0) x 107? at 2o level; (3) The LEP-I bound on the invisible Z-decay, 


T(Z — X) < 1.76 MeV, and the LEP-II upper bound on c(e*e^ > x9X9), which 
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is 5 x 107? pb for i > 1, as well as the lower mass bounds on the sparticles from the 
direct searches at LEP and the Tevatron; (4) The constraints from the direct search for the 
Higgs bosons at LEP-II, including the decay modes h > hi4hi,a1a4 — 4f, which limit all 
possible channels for the production of the Higgs bosons; (5) The constraints from B-physics 
observables like B > X,y, B, —^ pty, Bt — rtv, Y — yay, the ai mixing and the 
mass difference AM, and AM,; (6) The newest results for Higgs, top and stop results of the 
LHC. These constraints have been encoded in the package NMSSMTools [11]. In addition to 
the above experimental limits, we also consider the constraint from the stability of the Higgs 
potential, which requires that the physical vacuum of the Higgs potential with non-vanishing 


vevs of Higgs scalars should be lower than any local minima. 


130 mer me E S. rr —— 1.1 — 
128 L ^ with right-handed neutrino 4 = E. us 
x without right-haned neutrino 1.08 - * 
126 - * = s m : 
< e o 
m 1.06 - 
s1 E 
> > 
Q 12- I 
© E 104 - 
Ne 
a120 = 
E D use 
118 - 
" =. 
L n 
met 4 £ if 
dac * 5 E à 
L i a ^] 0.98 - 
112 [a ba doacaca d acaca dl acaca Lr saca T eger pe ques Au 1e 
500 600 700 800 900 1000 114 116 118 120 122 124 126 
A, (GeV) m, NMSSM«RHN (GeV) 


FIG. 8: The left panel shows the loop-corrected mass of the SM-like Higgs with or with- 
out the right-handed neutrino/sneutrino contribution, while the right panel shows the ratio 


NMSSM+RHN 


es xL NMSSM with my (mn MSSM) denoting the SM-like Higgs mass with 


(without) the right-handed neutrino/sneutrino contribution. 


The numerical results of our scan are shown in Fig. 8 in which we show the SM- 


like Higgs mass versus the tri-linear parameter A; in the left panel, and the ratio of 


pp CEPR / 77 NMSSM po PN 


versus mM in the right panel. Again we see that the contri- 
bution of the right-handed neutrino/sneutrino is sizable, which helps to push up the SM-like 
Higgs boson mass and thus makes the NMSSM more natural. 

Note that from Figs.6 and 7 we see that the correction to the Higgs mass can be positive 


or negative, depending on the parameter space (the most sensitive parameter is A). However, 
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under current experimental constraints the results in Fig.8 show that in the major part of the 
survived parameter space the correction is positive. The reason is that the parameter samples 
which give negative corrections are hard to survive the current experimental constraints 


(especially the Higgs mass lower bound given by LEP-II). 


V. SUMMARY 


In order to have massive neutrinos, the right-handed neutrino/sneutrino superfield must 
introduced in SUSY. In the framework of NMSSM such an extention will dynamically lead 
to a TeV-scale Majorana mass for the right-handed neutrino. Further, through the cou- 
plings SNN and SH,H,, the SM-like Higgs boson can naturally couple with such TeV-scale 
right-handed neutrino/sneutrino. As a result, the right-handed neutrino/sneutrino may 
significantly contribute to the Higgs boson mass. In this work we performed an explicit 
calculation and found that the Higgs boson mass can indeed be sizably altered by the right- 
handed neutrino/sneutrino. Such new contribution can help to push up the SM-like Higgs 


boson mass and thus make the NMSSM more natural. 
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Appendix 


Here we list the analytical renormalized formula for the elements of the Higgs mass 
matrix. Although they can be found in Ref. [12], we checked them and modified them 
according to our convention. Note that 6g denotes the tree-level 8 and the cx, sx, tx 
denote respectively cos X, sin X and tan X. 


The scalar 3 x 3 mass matrix Mà in the basis ^? = (Hu, Ha, S) is given by the entries 
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€CcgCg 
Mz, = 2Mwswd Mentes + tH, 88, (tgtg, + 2)] 
c2 932 M2. S2. c2 
t [MP + (M26 — M?,)c4,] + 4 (35) 
2 2472 
2 1,778 6 Beg 2 SBCB ry 42 2 242 A^ Mw sw $28 
M$., = IMwsw, tu tets, T tg, | = aa [M gr d (M7 == Mwy )cagl + e2 ? (36) 
2 2 2 
C2C V 2Mw swsec 
M, = Pe s gh, + try] + oO Ug chp Mau] 
dE E CUsCAg 
2AM; 2 —2 V2 X M3, 63 sac? 
_v2 WSwCgU PUES wewsetp (37) 
€ 8 €?v, 
2 
€C3C5 
M$,, = Mw swoh, a (2tetg, + 1) E: tg, tg] (38) 
82 2X2 M2, den 
+ > |M} + (MZtg? — Midàs + 1E, (39) 
2 2 
$8CB8C5 V2Mwsws2cg 
M$, = N LL + tH] + EUER - Mis] (40) 
sCAB s 
2AM i —24/2)A? M3 s3,,s2c 
V2 WwSwCgU ORE /2 Mir WEBS. (41) 
€ €3v, 


ty, Mwswsec3 
Viv, PUR, 


[2\? Misty s2a — 2kAe?v? — Mie sog] . (42) 


Mi. = KA,Us + Akt? EE 


[2M}+ Mwswsg — e(tn,tgso, + tn.c5,)] 


2 2 
My Sty S28 


m 2e*y2 


The entries Mp,, = Mp, (i, j = 1,2,3) of the pseudoscalar 3 x 3 mass matrix M2 in the 


i 


basis h? = (a, as, G)7 read 


20? Mi Si CA 
Mè, E o + Mg 5 My CA p; (43) 
Mw sw 828 cech 
M2 EE ce Me — —= [tg +t tot? 
Pi2 V2evscag H W Ag] xL Ay Ha" 3p] 
AM 
+ aan 2 M Siy S28 — 6ke^v?], (44) 
MÀ. = Mata + MAP 19262 ua scale 45 
P H= n —53 | ELE T EAE BR Hal ( ) 
tH Mwsw $8656; 
2 E s B 2 
MP, = -3A&KV, + v," Gua UU + tgstats,] 
M2, s352 AM2, s2 
Se 25 2p [Mz — My CA g] + STENT E Bap + 3ke^v?], (46) 
ev2cA g ety? 
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[1] 
[2] 


2 Mwsws2g 2 2 C8 C5, LAB 2 
—— —— D Mrg-tAg — Mi $2A^g] — lta, + tas teta, 


"^^^ 2V2ev.cag V2usCag 
AM 
+ Tah Mir Sr — 6Ke7v?], (47) 
M3, sin? AB 
Mj, = Ms tan? AB + aS —e’| 
€ 
[Lu,Cg-285 — in,Sg-a8s] - (48) 


2Mw swag 


J. Wess and B. Zumino, Nucl. Phys. B70 (1974) 39; M.F. Sohnius, Phys. Rep. 128 (1985) 39; 
S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150; N. Sakai, Z. Phys. C11 (1981) 
153. 

For a review, see, e.g., H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

J. E. Kim and H.P. Nilles, Phys. Lett. B138 (1984) 150. 

See, e.g., J. R. Ellis, et al., Phys. Rev. D39, 844 (1989); M. Drees, Int. J. Mod. Phys. A 
4, 3635 (1989); P. N. Pandita, Phys. Lett. B318, 338 (1993);Phys. Rev. D50, 571 (1994); 
S. F. King, P. L. White, Phys. Rev. D52, 4183 (1995); B. Ananthanarayan, P.N. Pandita, 
Phys. Lett. B353, 70 (1995); Phys. Lett. B371, 245 (1996); B. A. Dobrescu, K. T. Matchev, 
JHEP 0009, 031 (2000); V. Barger, et aL, Phys. Rev. D73,(2006) 115010; R. Dermisek, 
J. F. Gunion, Phys. Rev. Lett. 95, 041801 (2005); G. Hiller, Phys. Rev. D70, 034018 (2004); 
F. Domingo, U. Ellwanger, JHEP 0712, 090 (2007); Z. Heng, et al., Phys. Rev. D77, 095012 
(2008); R. N. Hodgkinson, A. Pilaftsis, Phys. Rev. D76, 015007 (2007); Phys. Rev. D78, 
075004 (2008); W. Wang, Z. Xiong, J. M. Yang, Phys. Lett. B680, 167 (2009); J. Cao, J. M. 
Yang, JHEP 0812, 006 (2008); Phys. Rev. D78, 115001 (2008); U. Ellwanger, C. Hugonie 
and A. M. Teixeira, Phys. Rept. 496, 1 (2010); J. Cao, Z. Heng and J. M. Yang, JHEP 1011, 
110 (2010); M. Maniatis, Int. J. Mod. Phys. A25 (2010) 3505; U. Ellwanger, Eur. Phys. J. 
C 71, 1782 (2011); J. Cao, et al., JHEP 1206, 145 (2012); JHEP 1203, 086 (2012); JHEP 
1210, 079 (2012); Phys. Lett. B 703, 462 (2011); arXiv:1303.2426 [hep-ph]. 

P. Fayet Nucl. Phys. B90 (1975) 104, Physics Letters Volume 69B, 4. 

D. G. Cerdeno, C. Muñoz and O. Seto, Phys. Rev. D 79, 023510 (2009). 

R. Kitano and K. y. Oda, Phys. Rev. D 61 (2000) 113001. 

D. G. Cerdeno and O. Seto, JCAP 0908, 032 (2009). 


17 


10 
11 
12 
13 
14 


J. Cao and J. M. Yang, Phys. Rev. D 71, 111701 (2005) [hep-ph/0412315]. 

U. Ellwanger et al., JHEP 0502, 066 (2005). 

K. Ender, T. Graf, M. Muhlleitner and H. Rzehak, Phys. Rev. D 85, 075024 (2012). 
D.J. Miller, R. Nevzorov, P.M. Zerwas, Nucl. Phys. B681, 3 (2004). 

K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G 37, 075021 (2010). 


18 


